Scaling of the quasiparticle spectrum for d-wave superconductors 
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In finite magnetic field H, the excitation spectrum of the low energy quasiparticles in a 2 dimen- 
sional d-wave superconductor exhibits a scaling with respect to H 1 ^ 2 . This property can be used 
to calculate scaling relations for various physical quantities at low temperature T. As examples, 
we make predictions for the scaling behavior of the finite magnetic field specific heat, quasiparticle 
magnetic susceptibility, optical conductivity tensor, and thermal conductivity tensor. These pre- 
dictions are compatible with existing experimental data. Most notably, the thermal Hall coefficient 
K xy measured by Krishana et al. in YBCO is found to scale as k is ~ T 2 F(aT/H 1/2 ) for T < 30A" 
in agreement with our predictions, where a is a constant and A is a scaling function. 
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The scientific community has been slowly coming to 
a consensus that the high T c superconductors have a 
d-wave order parameter One of the major differ- 
ences between these and conventional superconductors 
is that d-wave superconductors have gapless low energy 
excitations in certain directions in fc-space, whereas s- 
wave superconductors are gapped. As a result, the low 
temperature behavior in the d-wave case can be quite 
different from that of conventional superconducting ma- 
terials. Thus, in order to properly interpret experiments 
on these novel materials (and eventually develop a mi- 
croscopic theory), we must have a clear understanding of 
the physics associated with a d-wave order parameter. In 
this work, we attempt to elucidate some of this physics 
by deriving scaling relations obeyed by the quasiparticle 
energies and eigenfunctions. Using these relations, we 
are then able to deduce scaling properties of a number of 
important physical quantities. 

We begin our analysis with the Bogolubov || equa- 
tions, Titp = eip where ip = (u, v) is a Nambu 2-spinor 
whose components are the particle- like and hole- like part 
of the quasiparticle wavefunction respectively. Here, 

(h + V-E ¥ A \ 

H= { A* -h*-V + E F ) (1) 

with E F the Fermi energy, and h the kinetic part of the 
effective single particle Hamiltonian, and V(r) is the ef- 
fective disorder potential. As a model system we will 
choose to work with the effective one-particle Hamilto- 
nian h — ^ p „ A - ) with m the electron effective mass. 

2m 

Although this neglects any direct Hartree or exchange 
pieces of the interaction, such pieces are thought to be 
relatively unimportant except in renormalizing m and V. 
(In this paper we have set the charge of the electron e, 
the speed of light c, and Planck's constant fi all to unity.) 

In Eq. [j], A is the gap operator for spin singlet su- 
perconductivity defined as Ag(r) = J dr'A(r, r')g(r') for 
any g(r), where A(r,r') = — v(r — r'){ip (r)ip, (r')) with 
v the inter-electron interaction. If we rewrite A(r, r') in 
terms of center of mass coordinate R = ^iL. anc [ rela- 
tive coordinate x = r — r', then Fourier transform with 
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respect to x, we can write the gap function as A(R, k). 

In this work we consider a two dimensional d-wave 
superconductor. It is believed that this accurately rep- 
resents the high T c materials. We choose to consider 
a gap function with pure d xy symmetry rather than 
d x 2_ y 2 for notational simplicity. The final results for 
d x 2_ y 2 are identical. The gap function is written as 
A(R,k) = Ad Ty (R)k x k y j '(fc F ) 2 . Shifting back to the co- 
ordinates r and r', then integrating by parts, the gap 
operator can be re-expressed as 

s = iK.K.^W}} (2) 

where p x and p y are the components of the momentum 
operator, p F is the Fermi momentum, and here the brack- 
ets represent the symmetrization, {a, 6} = ^(ab + ba). 
The function Ad XH is the d-wave order parameter used 
in Ginzburg-Landau theory We can then consider 
calculating Ad XH in an inhomogeneous system by using 
a Ginzburg-Landau approach, then using A^ in Eq. [j] 
to find the quasiparticle spectrum. We note that this 
approach is not fully self-consistent in the sense that we 
will not use the derived quasiparticle states to then re- 
calculate the gap function. 

For a homogeneous system there are gapless nodes 
on the Fermi surface at the points p = (±p p ,0) and 
p = (0, ±p F ) where A vanishes. To study the low lying 
excitations near these points, we linearize the Hamilto- 
nian. As an example, we consider linearizing around the 
point p = (p F ,0). We write ip = e FX tp such that we 
can recast the Bogolubov equations as (Ho + H^ip = tip 
where 7io is the leading linearized term 

_{ Vv( Px -A x ) + V ±{ Py ,A d3sy (r)} \ 
n °-[±{ Pv ,A dxv *(T)} Vf (- Px -A x )-VJ [6) 

and Tii is the remaining piece 

where v F = p F /m is the Fermi velocity. 
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For sufficiently small energy excitations, Ho is much 
greater than Hi and it will be a reasonable approxima- 
tion to neglect Hi. To determine when this is a good ap- 
proximation, we consider the homogeneous case of Ad xy 
a real constant with A = 0, and V = 0. We then 
find that Ho is just the Dirac Hamiltonian for mass- 
less fermions in two dimensions, and thus has a coni- 
cal linear spectrum of quasiparticles e p = ±£(p) with 
£(p) = ^{v F p x ) 2 + (A dmy p y / p F ) 2 ■ Note that the coni- 
cal spectrum is highly anisotropic since v F 3> jp F . 
For excitations at temperature T, the typical momenta 
are p x ~ T/v F and p y ~ Tp F j Ad my - The largest term in 
Hi is then the term p 2 /(2m) which would be on order 
E F (T I Ad m ) 2 . For YBCO, photoemission spectroscopy 
H indicates that E F « 3,000X, and A dxy w A50K. 
In YBCO, we also note that the Fermi surface is not 
circular, but is somewhat flattened at the nodes (more 
square-like with rounded corners). This means we should 
really use an effective mass m y in the p 2 /(2m) term of 
H\ which lowers the energy scale of Hi by another fac- 
tor of perhaps two or three. Thus, we estimate that 
Hi /Ho ~ T/(2WK), so that the condition H > Hi 
may be well satisfied at temperatures as high as 30A'. 

We note that in Eq. [| a nonzero A x acts as a scalar 
potential for the quasiparticles. If we consider the case 
of a uniform superfluid velocity, we can choose a gauge 
where A<j is a real constant and A is proportional to 
the velocity (i.e., London gauge). In this case, A x in Eq. 
H acts as a scalar potential to yield a Doppler shifted 
spectrum e p = v F A x ± £ (p). Of course, were we to con- 
sider quasiparticles near the opposite (— p F , 0) node, we 
would have a spectrum e p = — v F A x ±£(p) . As pointed 
out by Volovik ||, in a magnetic field, screening currents 
exist with a typical velocity proportional to H 1 / 2 so that 
in a semiclassical approximation, the density of states at 
zero energy is proportional to H 1 / 2 . 

Let us now apply a magnetic field H perpendicular to 
the plane of the sample to create a vortex lattice such 
that the phase of the gap Ad twists a full 2ir as we 
go around each vortex. (Since the screening length is 
very long, we can assume H is homogeneous). The dis- 
tance between the vortices is proportional to the mag- 
netic length Ih ~ H~ 1/>2 . We now claim that at low T, 
to a very good approximation, the Hamiltonian Ho has 
a simple scaling form that we write (in a slight abuse 
of notation) as H$(r) = [H / ' H ]^H" c >(t[H / 'H ]%). In 
other words, if we can find the eigenvectors (r) and 
eigenenergies of the Hamiltonian Ho in field Ho, then 
the eigenenergies and eigenvectors in field H can be writ- 
ten as 

ff(r)=ffi(T[H/H ]t) (5) 
e%=[H/H ]h%°. (6) 

The first of these equations is the statement that the 
functional form of the eigenvector scales as the vortex lat- 
tice, whereas the second is just a reflection of the Hamil- 



tonian being linear in momentum. In order to demon- 
strate that these scaling properties hold, we consider each 
term in Ho individually. It easy to show that the vec- 
tor potential in a field H can be written in a scaling 
form A/j-(r) = [H/H ]i A Ho (r[H/ H ]^), and similarly, 
p must scale as the inverse of the characteristic length 
Ih so that p H = [H/Ho]^p H °- Thus we need only ex- 
amine V and Ad . 

We first consider the scaling of Ad xy ■ In order to have 
the desired scaling of Ho, we must have Ad xy H (r) = 
Ad xy Ho (r[H/Ho]^)- This is simply the statement that, 
like the wavefunction ip, the functional dependence of 
Ad xy on position scales with the vortex lattice. This 
assumption is clearly not correct near the vortex cores 
where there is some fixed width to the core that is rel- 
atively independent of field. However, away from the 
cores, for T <C T c , the magnitude of Ad x is fixed 
and only the phase varies. For |A<j | fixed, it can be 
shown that the function form of the phase that mini- 
mizes the Ginzburg-Landau free energy has the correct 
scaling properties. Thus, for T <C T c , we expect that 
Hamiltonian will scale to a very good approximation. 

We note that in general, currents in d-wave supercon- 
ductors can introduce some small amount of s-wave com- 
ponent (A s ) to the order parameter [Q. It can also be 
shown that the above scaling laws hold to linear order in 
Ag/Ad^y , so that for Aj > A s the scaling laws (Eqs. 
H and ||) are preserved 0. 

Finally, we turn to consider the disorder term V. For 
Gaussian delta-function correlated disorder such that 
(V) = 0, and (V(r)V(r')) = V S(r - r'), the disor- 
der does not define a length scale so that given a re- 
alization of disorder V(r), another configuration V'(r) = 
[H/H ]^V(r[H/Ho]^) is equally likely. In other words, 
the disorder term (in an ensemble average) has the proper 
scaling properties to preserve Eqns. ^| and |[ For more 
general disorder with a nonzero correlation length, these 
scaling laws will no longer hold precisely. We note, how- 
ever, that since there is no Anderson's theorem 0] for d- 
wave superconductors, the introduction of disorder will 
reduce the overall value of the gap and thereby reduce the 
maximum temperature at which the condition Ho 3> Hi 
is satisfied. 

Neglecting disorder once again, A^ x and A can both 
be considered to be periodic functions with the periodic- 
ity of the vortex lattice 0. Due to this periodicity, the 
eigenstates can be divided into Brillouin zones with one 
band of excitations per zone. The first zone should have a 
maximum k- vector of approximately |k max | w l H , with 
Ih the magnetic length. The number of different zones 
with momentum less than some k is roughly (k / k max ) 2 . 
The typical energy scale of an excitation of wavevector 
k is ky/ Ad xy v F /p F . Thus, the typical energy E n of the 
n th band is given roughly by E n — k max ^fnA dxy /m ~ 
^JnuJcAd^y where uj c — B /m is the cyclotron frequency. 
It is then convenient to define the dimensionful con- 
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stant a = ^/m/Aj^ such that E 2 w noT 2 B. Finally, 
it will be useful to define the dimensionless parameter 
x = a T/H 1 ' 2 which is roughly the number squared of 
bands that are considerably occupied at temperature T. 
For YBCO, a w .05Tesla 1/2 /^- 

Using the above described scaling laws we can extract 
a number of important statements about physical quan- 
tities. As a first example we examine the specific heat. 
We write the energy as 

n 

= [H/H Q ]^Y,^°f{<°[H/HS /T) (7) 

n 

where / is the Fermi function. The volume v of the sys- 
tem here scales as 1% ~ 1/ff so that v = v^[H^/H]. 
Thus, the energy density can be written as U/v = 
H 3 / 2 Fu(aT / H 1 / 2 ) where F v is some scaling function 
that we can write down in terms of eigenenergies, but can 
not completely evaluate without fully diagonalizing 7io- 
Here we have used the fact that the sum in Eq. ^ is only 
a function of the dimensionless quantity x = aT/H 1 / 2 . 
Note that here and elsewhere in this paper, we will as- 
sume that we are sufficiently far below T c that the mag- 
nitude of the gap does not change much with T. Differ- 
entiating to obtain the electronic specific heat per unit 
volume we find that 

C v =TH^F c (aT/H%) (8) 

where Fc is again some unknown scaling function. Wc 
note that this scaling form does not include contributions 
to the specific heat from normal state electrons in the vor- 
tex cores. These contributions, however, are thought to 
be small ||. As discussed above, in a semiclassical ap- 
proximation, the density of states at zero energy is pro- 
portional to if 1 ' 2 so that the specific heat is proportional 
to Tif 1 / 2 at low temperatures as predicted by Volovik || . 
Thus the function Fc should be a constant at small ar- 
gument. Deviations from this constant should first occur 
when we start to fill more than one band. This happens 
at x w 1, or at (T/ff 1 / 2 ) w 20 if/Tesla 1/2 . 

Similarly, starting with the free energy F = U — TS, 

where S = J2 n (/« ln /« + 0- _ fn) m (l ~ fn)) and f n = 
f(e n /T) we can write the free energy per unit vol- 
ume in the scaling form Fjv = if 3 / 2 F F (T/ff 1 / 2 ). We 
then conclude that the quasiparticle magnetic suscepti- 
bility per unit volume scales as x — d 2 (F/v)/dM 2 = 
hF x (aT/H^) with F x an unknown scaling function. 
Since there is no crossing of states through the Fermi 
level as we change magnetic field, we do not predict any 
de Haas- van Alphen oscillations. However, there may be 
oscillatory contributions to the susceptibility from the 
condensed fraction and the normal vortex cores that we 
do not consider here 

We now turn to consider electrical and thermal trans- 
port properties. We first define [O the charge velocity 



operator v' 1 ) = i[H, ra z ] and the thermal velocity opera- 
tor v^ 2 ) = i[H, {H,r}]. Operating on a state ifi — e lkFX tp 
near the node at (p F ,0), we find 

v (i) e ik F x^ = e ifcFX^ Q) { azj r}]$ + smaller terms (9) 
v (2) e *fc F *^ = e «F«[7i 0| {H ,r}}^ + smaller terms (10) 

where the smaller terms are typically smaller by order 
"Hi /"Ho- It is then easy to see from this form that the 
operator v^ 2 ) scales as if 1 / 2 whereas the operator 
scales as if . 

We now use the Kubo formula to write the generalized 
response function at frequency to by pd[ | 

Lab _ T ^ (n|v< a) |m)(m|vfV)/( £K /r, £m /T) 
ij v ^ (e„ - e m - w - i0+)(e n - e m + i0+) 

where / is the thermal occupation factor, v is the volume 
of the system, the indices i and j take the values x and y, 
and the indices a and b take the values 1 and 2 (for charge 
and heat transport respectively) . Noting that the volume 
of the system scales as if -1 and the energies all scale as 
if 1 / 2 we immediately obtain the two parameter scaling 
law - T a+b - 1 F^ b (aT/H 1 ^, aw/H 1 ' 2 ), where F£ b 
is again some scaling function that we will not evaluate. 
The real part of the optical conductivity tensor is defined 
as Re [cry ] = ^RefL 11 ] which immediately yields a two 
parameter scaling law for the optical conductivity. 

Re[ofc] ~ F^iaT/H 1 ' 2 , auj/H 1 ' 2 ). (11) 

It should be noted that in this Kubo formula calculation 
the response of the superfluid fraction has been neglected. 
This then does not include, for example, the response of 
the system due to the motion of vortices |l2f . 

We now turn our attention to the DC (w = 0) ther- 
mal conductivity tensor k, defined as the matrix that 
relates the heat current j 9 to the temperature gradient 
via j q — kVT '. We note that experimentally, a large 
part of the diagonal components of this tensor are due 
to phonon transport of heat. However, the Hall (off di- 
agonal) component of this tensor should be completely 
electronic in origin [0. Note, that when calculating re, 
one must usually take into account the effect of the ther- 
moelectric coefficient L 12 . However, here we can neglect 
that contribution, since there is never any voltage in the 
superconducting state. Thus, we have rey = t^L 22 , and 
we obtain the naive scaling law re y - TF 22 (aT / H 1 ' 2 ). 
Although this is indeed the correct scaling form for the 
(electronic part of the) diagonal component of the tensor, 
it is not correct for the Hall component. It can in fact 
be shown that the scaling function F 22 here is precisely 
zero due to the particlc-holc symmetry inherent in the 
linearized Hamiltonian 7f-o- This result is very easy to 
understand. Imposing a heat source on one side of the 
system excites many particles and holes. Both particles 
and holes diffuse in the direction of the heat sink. In a 
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FIG. 1. Thermal Hall data from reference (l|. (a) Ther- 
mal Hall transport coefficient K xy plotted against external 
magnetic field H at temperatures (from bottom to top) 
20K, 22.5JC, 25K, 27.5K, and 30K. For technical reasons data 
has not yet been taken below 20K. (b) Collapse of these 
five curves according to the scaling law shown in Eq. 
Note that the characteristic scale of H 1 ^ 2 /T is approximately 
.05 Tes\& 1/2 /K. 



Al Hk \ - \m H ){m H \Hi\n H ) lTT/TT ,i/2 x , h oX 

m " m 

Note that both of these corrections scale as H 1 / 2 . In- 
cluding these first order corrections into the Kubo for- 
mula and expanding we find these correction terms give 
a leading contribution to the thermal Hall conductivity 
that scales as 

Kxy ~T 2 F Rxy {aT/H l l 2 ) (12) 

with F Kxy again some scaling function. As shown in Fig. 
1, experimental results of reference [fLs)] do indeed show 
this scaling form at temperatures below 30K (For tech- 
nical reasons data has not yet been taken at tempera- 
tures below 20K). As discussed above, the characteristic 
scale for features in the function F^ v (i.e., where the 
curve becomes nonlinear) should be seen at x w 1 or 
H x l 2 jT w .05 Tesla 1/2 /^, which is in good agreement 
with experiment. 

In conclusion, we have found that the scaling proper- 
ties of the quasiparticle spectrum in d-wave superconduc- 
tors provides a very general and powerful tool for ana- 
lyzing various physical quantities. 
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L. S. Levitov. The authors are greatly indebted to N. P. 
Ong and K. Krishana for sending us their experimental 
data prior to publication. This work was supported by 
NSF Grant Number DMR-9523361. 



magnetic field, the particles curve one way and the holes 
curve the other way. Thus, when there is particle-hole 
symmetry, there is no net Hall transport of heat. 

The proof that the linearized Hamiltonian yields n xy — 
is a little bit involved. The particle-hole symmetry is 
expressed mathematically by saying that given an eigen- 
pair e n ,ip n satisfying TCoipn = enVVi, it can be shown 
that there exists another eigenvector ip' n = CF y ip^ with the 
same eigenvalue also satisfying TCoi^'n = £ n ^|, where a y 
is the usual Pauli spin matrix. Using this symmetry and 
adding up the contributions to n xy from all four nodes 
on the Fermi surface, it can be shown that n xy vanishes. 
It should be noted that n xy remains zero even when we 
include the smaller terms in from Eq. [l0| We also note 
that this mechanism does not force <x x „ = 0. Further 
details of this calculation will be published elsewhere J7[ . 

In order to find a nonzero n xy , we must break the 
particle-hole symmetry by including the contributions 
from Tii. As mentioned above, at low enough T, we have 
Tii <C tCo so that this term can be treated perturbatively. 
Inclusion of this term then shifts the eigenenergies via 
e„ — > e n + Se n and the eigenstates |n) —* |n) + 8\n). To 
lowest order these shifts are given by the usual expres- 
sions from perturbation theory which then obtain the 
following scaling forms 

Se» = (n H \H 1 \n H ) = [H/H } 1 / 2 5e^ 
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